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CONVERGENCE TO EQUILIBRIUM OF A BODY MOVING IN A KINETIC SEA 


XUWEN CHEN AND WALTER A. STRAUSS 


Abstract. We consider a continuum of particles that are acted upon by an external force G(t, x) and 
that collide with a rigid body. The body itself is subject to a constant force E as well as to the collective 
force of interaction with the particles. We assume that the particles that collide with the body reflect 
probabilistically with some probablility distribution K(v,u). Under certain conditions on G(t, x) and 
K(v,u ), we identify an equilibrium velocity Voo of the body and we prove that this equilibrium is 
asymptotically stable. 


1. Introduction 

In previous work of this type the particles have been subject to no force (G(f, x) = 0) and therefore 
travel with constant velocities between their collisions with the rigid body. In this paper we initiate 
an investigation of how some forces on the particles may affect the motion of the body. What is the 
asymptotic effect of such particle forces on the body? Our ultimate goal for the future is to treat a 
plasma in which the particles are charged and thus are subject to electromagnetic forces, which may be 
either external or created by the particles themselves. In this paper the forces are external. 

Our problem has a free boundary, the moving location of the body. The other unknowns are the 
configuration and motion of the particles. The particles may collide with the body diffusely. Boundary 
interactions in kinetic theory are very poorly understood, even when the boundaries are fixed. Free 
boundaries are even more difficult. Our most important assumption is that the whole system, consisting 
of the body and the particles, starts out rather close to an equilibrium state. 

We consider classical particles that are extremely numerous and subject to an outside force G, which 
in this introductory study we assume to be small and decaying. While one could consider modeling the 
particles as a fluid, we instead model them as a continuum with a phase-space density / as in kinetic 
theory mi mens- Furthermore, the interaction of the particles with the body at its boundary can 
be quite complicated in typical physical scenarios. For instance, the boundary may be so rough that a 
particle may reflect from it in an essentially random way. There could even be some kind of physical or 
chemical reaction between the particle and the molecules of the body. Therefore we model the collisions 
by a probability distribution K. This is sometimes called diffuse reflection. Altogether, the prescribed 
data consists of the initial density /o, the force G, the collision distribution K 1 and a constant force E 
on the body. 

The present paper is a sequel to [3 [El and is also highly motivated by a series of remarkable papers 
ID 0 0. Most of the earlier papers, including ID 0 El El llOj, were devoted to proving more detailed 
behavior such as determining the precise rate of approach to equilibrium. In particular, the earliest 
papers 013 assumed specular reflection at the free boundary. On the other hand, considering a perfect 
gas, the authors of chose a Maxwellian-like distribution for the reflected particles. More generally, 
for the reasons mentioned above we considered in 0H3 as well as in this paper a general probabilistic 
law of reflection. In 0113, where G = 0, we found a condition on K and /o that is almost necessary 
and sufficient that the approach to equilibrium is reversed; that is, the approach begins from the right 
side, say, and ultimately approaches from the left. The present paper does not deal with such detailed 
behavior, focusing instead on whether it approaches equilibrium at all. Some numerical investigations, 
mostly dealing with a diffuse reflection which ejects Maxwellians, have appeared in ElDSUIZIQSl- Some 
other related investigations that deal with the speed of approach to equilibrium in the specular reflection 
case are m 0 0 E3. See also the survey in [3]. In all of the preceding references, the particles are 
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identical, the body is initially moving near the equilibrium velocity V» and the motion of the body is 
one-dimensional, even though the particles move in three (or d) dimensions. The current paper is the 
first one for which the particles are themselves subject to a force G/0. 

We now describe our problem in mathematical terms. For simplicity we take the body to be a vertical 
disk D{t) centered at the point {X(t), 0,0) and traveling with velocity (V(f),0,0). We assume there is 
a constant horizontal force Egl acting on the body, in addition to the horizontal frictional force F(t) 
due to all the colliding particles at time t. Thus 


dX 

dt 


m = e-f( t). 


The particle distribution /(f,x, v) satisfies the Vlasov equation 


dtf + v • V x / + G(f, x) • V v / = 0 (1.1) 

in R 3 \D(f). We assume the initial velocity /(0,x, v) = / 0 (v) depends only on v and is even. At each 
collision with the body we assume the diffuse law of reflection 

/+(i,x,v) = f K(v! - V(t);ui - V(t)) f-(t,-x.,u 1 ,v±) dm (1.2) 

J (u 1 -V(t))(t> 1 -V(t))<0 

for x eD(t), where the pre/postcollision distributions f± are defined as 


f±(t,*,v)= lim f(t ±£,x±ev,v). 

£—> 0 + 

In (11.211 we regard u as the velocity of a particle coming into collision and v as an ejected velocity. We 
assume that the collision kernel K conserves mass during collisions. The form of OI implies that there 
is no energy exchange between the disk and the particles except in the horizontal direction. This is 
natural because the disk is fixed to move only along the z-axis. More detailed assumptions on G, /o and 
I\ are provided in Section [2j We note that the F(t) notation for the frictional force does not reflect the 
complicated nature of the interaction. In fact, due to the recollisions between the body and the particles, 
we will see in Lemma HOI that the frictional force at time t depends on the behavior of the whole system, 
including both the disk and the particles, at all previous times. 

It is far from obvious what the ultimate velocity V, x will be. We prove that Vx is entirely determined 
by imagining a situation in which both G = 0 and there are no precollisions; that is, the collisions at a 
time t come from particles moving at constant velocities that have not previously collided with the body. 
The limiting velocity is determined as follows. We define a “fictitious” force Fqo for which there is 
no external field G at all and the particles do not hit the boundary before time f, namely, 

F 0 o{V)= f dS x [ L(v! - V(i)) / 0 (v) dv, (1.3) 

JD(t) Jr 3 

where L is defined in (12.811 . It is a strictly monotone function of V. We then define Vx, by the equation 


FooiVx) = E. 

Our main result is stated informally as follows. If the initial velocity V(0) is sufficiently close to Voo, 
then there exists a solution (V(f), f(t, x, v)) of our problem such that V(t) —> Voo as t —> oo. The same 
asymptotic statement is valid for any solution of the problem. This result is stated with precision and 
further details in Section [2] The present result is the first one for which the particles are themselves 
subject to a force G ^ 0, so that between collisions with the body the individual particles travel in 
curves rather than straight lines. As a consequence, we have limited our investigation to proving the 
convergence (at some rate) to the equilibrium velocity Vx,. 

We state the basic iteration scheme to obtain the main theorem in (12.181) . The core of the proof is a 
series of technical estimates. We emphasize that adding an external force on the continuum of particles, 
even considering only the rather ideal situation where the external force G(f, x) acting on them is small 
and decaying, makes the problem much more difficult to deal with. Because we are interested in the 
long-time behavior of the body, even a small force G could have significant effects for large time which 
we explain below. 

Relating the aforementioned Vx>, which comes from an imaginary case in which there are no collisions 
and G = 0, to our system (for which there are collisions and G ^ 0) is the most delicate and novel part 
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of the analysis. In Section [3] we estimate the force H(t) on the body at a large time t due to G in the 
fictitious situation that there have been no precollisions. H(t) is given by an integral over the various 
velocities v of the incoming particles and we prove that H(t) is integrable over 0 < t < oo. We estimate 
the integral in five pieces, using an intermediate time 0 < T < t. One piece H i is small (for large t ) 
because v — is small. The second piece H 21 is proven small by using the decay of G in time. The 
third piece H 221 makes use of the decay of G in both space and time. The fourth piece H 222 makes use 
of the decay of G in space. The fifth piece H 223 particularly makes use of the decay of fo as well as 
the decay of G in time. There is a range of possible choices for the parameters p, m, q, summarized in 
Theorem 13.21 

Furthermore, because the body is moving and the particles travel in curves, it is evidently more 
difficult to track the collisions between the particles and the body, in contrast to the G = 0 case for 
which the particles travel in straight lines. A given particle can collide with the body many times or even 
an infinite number of times. This is partially discussed in Lemma 12.21 and at the beginning of Section 
[4] Nevertheless, the effect of the precollisions can be estimated by the most recent precollision. This 
delicate part of the proof is given in Lemma POl Finally, putting the different pieces together is a matter 
of dealing with several competing small quantities, as we discuss in Section [5j where we combine the 
previous estimates and prove the required convergence. 

1.1. Acknowledgements. The authors would like to thank Kazuo Aoki once again for introducing 
them to this general subject and Yan Guo for discussions about multiple collisions. This research was 
supported in part by NSF grant DMS-1464869. 


2. Assumptions and preliminaries 


We write the coordinates as x = (xi,a;j_) = (xi, X 2 , £ 3 ), v = (vi,i>_l) = (ui, V 2 , U 3 ), G = (Gi,G_l) = 
(G 1, G2, G3). We begin by listing the assumptions precisely. We assume that the force G has the form 

G(t, x) = (Gi(f, £ 1 ), G±{t, £_l)) (2.1) 


and that it decays like 

\Gi(t,xi)\ < and l G J-(M±)| < (2.2) 

where cg is a constant, q > 2 and m > 0. The structural assumption (EH) is crucial to our proof although 
we hope to weaken it in future work. In the main theorem below we obtain a decay rate that does not 
depend on the dimension because we refrain from making an even stronger assumption. In fact, even 
with our decay assumption E3 on Gj_, some of the particles might collide with the disk many times 
and might not even escape the disk at all. 

We assume that the initial velocity distribution /o is a nonnegative even G 1 function ^ 0 that has the 
form 

/o(v) = ao(ui) &o(u_l) (2.3) 

and that it decays like 

|V ’ / " WI 4 („)'■+>,/■+' <2 ' 4) 

where l\ > q + 1 and I 2 > 1- (We remark that it is sufficient to assume the slightly weaker condition 



(vi) q di>i < 00 , 



|V&o(?Ll)Mu± < 00 


although we will not bother to explicitly use this condition.) 

Define 7 = |V(0) — I4o|. This quantity is used throughout the paper and will be assumed to be very 
small. We assume that K(v\,u{) is a continuous nonnegative function, G 1 for u\ 7 ^ 0,'iq 7 ^ 0, that 
is even in both variables separately (K(vi,ui) = I\(—Vi,ui) = K(vi,— rti)) and that is bounded for 
v\ £ R, |'iii| < 3y. We impose the mass conservation condition 

v 1 K{v ll ui)dvi = \ui\, (2.5) 
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the power law condition 


pOO 

sup / ViK(vi, u\)dvi < C\u\\ p 

\<~i+cn J 0 


l“i I<7 +cg J 0 

with 0 < p < 2 , and the integrability condition 

K(vi,z- y - Voo) ( z)~ h < M(z) for |ui| < 37 , \y\ < 37 , \z\ < 00 
where M £ L 1 (R). We also define 

L{u\) = u\ + / vfK(vi,Ui)dvi and L(u\) = sgn(«i)L(ui). 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


Because K is even, L(ui), as defined, is also even. In addition, we assume that the even function L{u\) 
is decreasing for u\ < 0. Then we define the collision operator as /C t (/_) = /+ as in (11.21) . 

To summarize, the assumptions listed above have five exponents p, to, q, l\, I 2 that satisfy 0 < p < 
2, q > 2, to > 0, l\ > q + 1, I 2 > 1. We define another parameter a by 

- = p = min(TO, q — 1). (2.9) 

a p+1 p 

We also have the two parameters 7 and cg in (12.21) that will be chosen sufficiently small. Our main 
theorem is as follows. 


Theorem 2.1. Let p > 1 + There is a constant A such that if 7 = |V(0) — I4o| and cg are 
sufficiently small, then there exists a solution (V(t), f(t, x, v)) of our problem for which V £ C ll ([0, 00 )), 
/ € L°°([0, 00 ) x R 3 x R 3 ) and 

\V(t)-V x \<ye- b ° t + -^^ ( 2 . 10 ) 

where bo = miny € ry oo _ 3 7 j y oa+37 ] Fq 0 (V). The pair of functions /±(t,x, v) are a.e. defined explicitly in 
terms ofV(t) and /o(v). Furthermore, any solution of the problem (in the sense stated above) satisfies 

(I2T01) . 

We now mention three examples of collision kernels for which the reflected velocity distribution is 
Gaussian for each incoming particle. 

Example 1. Let 

= 2)3e~ Pv ' |ui|. 

As shown in jH, IQ , K{v\,u\) satisfies the assumptions with p = 1. In this case, we require q > 3 and 
to > 2 in 12.21) . and l± > 4 in |R^| ). 

Example 2. We now choose 

v 2 

K{vi,u\) = 2e 

It is proved in GO CGI] that this kernel checks our assumptions with p = |. In this case, we need q > 8/3 


and to > 5/3 in K2.2 1) . and l± > 11/2 in \2.4\ ). In this example an incoming particle with almost the same 
velocity as the body is likely to be reflected with almost the same velocity, while an incoming particle with 
a very different velocity is reflected according to a very wide Gaussian around V(t). 

Example 3. We now scale the last example by setting 

K(vi,ui) = C \ Ul f e~ v ^ u ^~ X 

where /3 £ [—1,3) and C is chosen such that mass is conserved, that is, 12.50 is satisfied. The most 
important feature of this kernel is that, as /3 runs through [—1,3), p runs through (0,2]. At the endpoints 
of this interval we have the following behavior. If p = 2, then q > 5/2 and to > 3/2 in S2.20 . l\ > 7/2 
in If p —>■ 0, then q, m — 00 in \2.20 . l± —>■ 00 in j2-4\ l- 


We proceed by pointing out some basic properties of our model. The total mass M = f dxj dvf(t, x, v) 
of the particles is an invariant. Indeed, it is obviously invariant under the flow (ED while it is also pre¬ 
served under collisions due to ED , as proven in [9} Lemma 2.3]. The total horizontal force at time t 
due to the collisions is 

dv L(v 1 -V(t)) f — (t, x, v) 


F(t) = / dS x 

JD{t) 


/ R3 


( 2 . 11 ) 
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as shown in [5] Lemma 2.2], It is the sum of the force due to particles on the right (iq < V(t)) and the 
left (vi > V(t)). 

Our basic technique is to look for possible velocities in the class W, where we define W £ W if and 
only if W(t) is Lipschitz, W(0) = Vo, and 

AjP +1 


\W(t)-V 00 \^'ye~ b ° t + 


( 2 . 12 ) 


(1 + t)° 

where the constant A will be specified later in Lemma 15.21 and the power cr is given by m • Among 
other conditions to be stated later, we require A > 1 and 7 << 1 so that A'y < 1 . For W £ W we define 
X(t) as the primitive of W(t): 

^ = W ; X(0) = 0. 

Until the last section of this paper we will fix any possible velocity W £ W. 

We introduce the characteristics of the PDE O) without collisions as follows. Given (f,x,v), define 
(x(s; t, x, v), v(s; t, x, v)) to be the solution to the characteristic equations 

dx(s;t,x, v) 


ds 

dv(s ; t, x, v) 
ds 


= v(s; t, x, v) 

= G(s, x(s; t, x, v)) 


(2.13) 


with final condition (x(f; t, x, v), v(i; t, x, v)) = (x, v). In particular, for G(f, x) of the separated form 
m, the horizontal and vertical components are decoupled: 


(£i(r;f,x,v),di(r;f,x,v)) = (xi{r-,t,x,vx),vx(r\t,x,v\)) 

(ij_(r;t,x, v),vj_(r;t.,x, v)) = (xj_(r; t, x±, v±), v±(r; t, x±, v±)). 

For any function <f> on phase space, we define the solution operator St r as 

(S tr . 0 )(x, v) = 0 (x(r; t, x, v), v(r; t, x, v)). 

Then of course /nb = Stofo solves the Vlasov equation 

(ft + vV x + G(t,x)-V v )/jvB = 0 

fNB (0, X, v) = / 0 (v) 


(2.14) 


(2.15) 


in K 3 with no boundary conditions. 

The fictitious force Fqq, previously defined in m, is a strictly monotone function of V, due to L 
being an even function that is decreasing for iti < 0. For completeness we state this as a lemma. As a 
consequence, V 30 will be uniquely defined by Foo(U 00 ) = E. 

Lemma 2.1. Suppose /o(v) ^ 0 is even, continuous and ^0. If L € C 1 is an even function with 
L( 0) = 0, and L'(u\) < 0 for u\ £ (— 00 ,0), then Foo(U) is an increasing odd C 1 function of V. 

Proof. A simple change of variable yields 

Foo(V) = off L(v 1 -V)f 0 (v)dv- [ L(vx-V)f 0 {v)dv) 

\Jv Jv i^V / 

= C ( f L(v 1 - V)fo(v)dv - [ L(v 1 + V)fo(v)dv\ 

\Jv i<V Jv i<-V / 

where C is a fixed constant. Thus Fqq is odd and its derivative is 

^oo(^) = ~C ( f L'(v\ - V)fa(y)dv + [ L\v x + V)f 0 (v)dv) > 0. 

\J Vl <V Jv x <-V / 


□ 


Given a time t, we will have to estimate the frictional force due to the collisions of the particles. To 
that end, we define another “fictitious” force as follows (though not as ’’fictitious” as Too). The fictitious 
force Fo(t) for which no particle hits the boundary before time t is 

F 0 (t) = F 0 (t,W(t)) = [ dS x [ L(vi - W(t)) f NB {t,x, v) dv 
J D(t) JR d 


(2.16) 
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where /jvs(i, x, v) solves the kinetic equation (12.151) with no boundary condition and is subject to the 
initial condition /atr( 0, x,v) = /o(v). As before, the fictitious force F 00 is defined so that there is no 
external field G at all and the particles do not hit the boundary before time t, namely, 

F 00 = F 00 ( W{t )) = [ dS x [ L(v i - W(t)) / 0 (v) dv. (2.17) 

JD(t) JR d 

The force due to the precollisions on the body moving at velocity W(t) then can be expressed as 


Rw{t ) 


F(t,W(t))-F 0 (W(t)) 





L(v 1 - W{t)) \f NB (t, x, v) - /_(i,x, v)] dvdS x 


+ [ [ L{vi - W{t)) [/_(t,x, v) - f NB (t,x, v)] dvdS x 

JD(i) Jv 1 ^W{t) 

R\v(t) + 


comprised of the forces on the right and left sides. 

Given any W G W, we define the function Vw(t) by the iteration scheme 

^ = E y F00( ^) ( V °o - Vw) - Rw (t) + F 00 (W) - F 0 (t, W), V w (0) = V 0 . (2.18) 

Note that a fixed point (that is, Vw = W = V) would satisfy 

^ = E y F ™y ] (Voo ~v ) - R v (t) + Foo(V ) - F 0 (t, V) = E - Rv(t) - Fo(t , y) 

and dX/dt = V. Thus any fixed point solves our problem. We will estimate the term Fqo(W) — Fo(t, W) 
in Section [3] and estimate the term Rw{t) in Section [H 

A key difficulty in our problem is that a typical particle may collide with the body many, or even 
infinitely many, times. The following lemma, written for convenience in the ID case, illustrates this 
difficulty. It explains how one can represent the solution by an expansion in terms of a finite number k 
of previous collisions. This expansion might never reach the initial datum if there are infinitely many 
collisions; it could happen in many ways. 


Lemma 2.2 (Integral Representation). Let ( x,v ) Z(t) where Z(t ) is the set of measure zero defined 
in the beginning of Section^ For brevity we write Bt r = St r ICt and J(s,x,v) = — G(s, x)/q(u). Given 
t > 0, we denote the particle collision times as t > t\ > t 2 > .... Then for arbitrary k ^ 2 we can 
represent the solution to (ED subject to the collision law \ 1 . 2 \ ) by 

m = r 5 tjS J( S )d S | (2.19) 

1 ( dmax(tj + i,0) J 

+ j S ts J(s)ds. 

J max(ii ,0) 

Proof. It is important to note that the actual number of precollisions could be larger than k in formula 
(12.191) . Thus if there are an infinite number of collisions, /o would never show up in (12.191) . Note also 
that each tj depends on t, x , v as well as on the pre/postcollision velocities at all L for 1 < * < j. Because 
the lemma is similar to j!2l Lemma 24], we only sketch the proof here. It is instructive to first write out 
some of the terms in formula (12.191) in total detail including all the variables. For example, we have 


B ttl g(x,v) = S t , tl ICg(x,v) = J dui K(v{t ^; t, x, v) - L(ti),wi - L(fi)) g(X(t i),ui) 


Bt tl B tl t 2 g{x,v) = J dux K(v(tf;t,x,v) - L(fi),ui - V(ti)) 
x J du 2 A"(z)(t^;ti, A(G),mi) - V(t 2 ),u 2 -V{t. 2 ))g(X(t 2 ),U 2 ) 
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where tf means lim. t +, that is, postcollision. In fact, exhibiting all the variables, we have 

B Ul ■ ■ .B tj _ ltj g(x,v) = J dui K(v(tf-,t,x,v) - V(t{),ui -V{ti)) x ... 

X I j tj— ii X(tj~ i), Uj— i) — V(tj),Uj — V(tj )) g(X(tj),Uj). 

We now begin the formal proof with the standard Duhamel formula 

f(t) = X{o<t 1 }Stt 1 >Ct 1 f-(ti)+ f S ts J{s)ds. 

J max(ti ,0) 

We then write the next Duhamel formula 

+ [ S tlS J(s)ds. 

J max(t2,0) 

Combining the two formulas, we obtain 

f(t) = X{o<t 2 }^tt 1 Bt 1 t 2 f-{t2) + £>«i f St lS J(s)ds+ f St s J(s)ds. 

«/max(t 2 ) 0 ) ^max(ti,0) 

This is the case k = 2. Iterating in the same manner by induction, that is, replacing f—(tj) with its 
expression at the next level in order to get /_(tj+i), we obtain formula (12.191) . □ 

3. Effect of the particle forces on the body 
The difference between the two fictitious forces (12.161) and (12.171) is 

H(t) = H = F 00 (t) — F 0 (t) = [ dS x [ L(vi-V(t))h{t,x.,v)dv 


1 D(t ) 


lR d 


where we denote 

h(t , x, v) = f NB (t, x, v) - / 0 (v) . 

In this expression there is no collision with the body. The function h solves 

(flt+v- V x + G(f,x) • V v )/i = —G(f, x) • V v /o (3.1) 

h( 0 ,x, v) = 0 

or written in Duhamel form, 

/i(f,x,v)= / S ts J(s)ds 

Jo 

where we also denote 

J(s) = J(s, x, v) = —G(s,x) • V v /o- 
For any T £ [0, t], we could also write 

h(t, x, v) = h(T, x(T ; t, x, v), v(T ; t, x, v)) - [ G(s, x(s; t, x, v)) • V v /o(v(s; t, x, v)) ds. 

Jt 

We claim that H (t) is small and integrable in time so long as G is small. To prove such a claim, we first 
need the following lemma. The constants denoted as C occurring in the estimates throughout this paper 
are independent of 7 , cg, t, and A as well as the solution. The constant A will be chosen sufficiently large 
relative to one of the constants C. 

Lemma 3.1. There is a constant C > 0 independent of cg, 7 , t, x, and v, such that 

'^ X ’ V) I < cgC 1 w(*i(s-,t,x,v)) mds {vi yi+i {v± )‘*+ 1 - CG ’ c ' (3 - 2) 


Proof. We have 


S ts J{s) = —G(s, x(s; t, x, v)) • (V v /o) (v(s; t, x, v)). 


So by assumption, 


|S ts J(s)| < c G C- 


(s) q (xi (s;t,x, v)) (i) 1 (s; t, x, v))* 1+1 (v± (s;t, x, v)) 


*2 + 1 ' 
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Hence 


|ft(t,x,v)|< [ \S ts J{s)\ds 
Jo 


(3-3) 


< c G C 

Finally notice that 


o { s ) q (*i( s ; ^ x > v )) (t>i(s; t, x, v ))* 1+1 (nj_(s; t, x, v )) i2+1 


ds. 


|^i(s; t, x, v) —Dil < j \G(p,Z(p;t,x,v))\dp ^2 J -^dp if c G C 
\v±(s;t,x,v) - u_l| < j \G(p,±(p;t,x,v))\dp ^ 2 J-^dp ^ c G C. 


(3-4) 


So we can replace vi(s; t, x, v) and hx(s; t, x, v) in estimate (El by v\ and v± at the price of a different 
constant C. Thereby we obtain estimate El- □ 

The next two lemmas provide rather delicate estimates of H(t). 

Lemma 3.2. We select T = t a with 0 < a < 1 to be determined later and assume 0 < p < 2 and the 
following (somewhat redundant) decay conditions: 


(1 — a) (p + 1 ) > 1 

(3.5) 

l\ > 2 and a {q — 1 ) > 1 

(3-6) 

am > 1, q > 2 

(3.7) 

l\ > q + 1 and a (q — 1 ) > 1 . 

(3.8) 

Then the two fictitious forces Fq and Fqo are near each other in the sense 

that there is C independent of 


t such that 


Proof. Recall 


\H(t)\ < c G C 


(t)(l-“)( P+1 ) ( t) am /A“(9-t) 


(3.9) 


H(t) = / dS x / L(vi — W(t))h(t,xi, v)dv. 

J D(t) J R d 

First we consider t < 1. Note from El that 

L(v r - W(i)) < C (|«i - W(t)\ p + \ Vl - W(t )| 2 ) < C( Vl ) 2 

because p <2 and 

|wi -W(t) | < |ui - Hoo| + |Roo - W(f)| < l^i| + |Fool + (7 + A l) < C'(ni) 

for t < 1. In the previous estimates we have used the assumption that sC 1 and 7 is small. By Lemma 
13.11 we therefore have 

\H(t)\< [ dS x [(v 1 ) 2 Cc G (vi)~ h ~ 1 (vj_y l2 ~ 1 dv < Cc G 

J D(t) J 

because l\ > 2 and I 2 > 1 . 

In the rest of the proof, we assume t > 1. We decompose H(t) = H\{t) + H 2 (t), where 

Hi (t) = f dS x ( L{v\ — W(t))h(t, x, v)dv 

J D(t) JS 


H 2 (t) = / dS x L(vi - W(t))h(t,x, v)dv 

JD(t) Js c 


S = i v GM d : \ Vl - Kol < ^ } . 


and 


The constant b will be specified later and we choose T = t a . 
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The estimate for Hi is fairly simple. In fact, we have by Lemma T3.II that 

1 


\Hi(t)\ < c G C 

< c G C 


dS x 


D(t) 


'\vx-VvoKHf 


dv\L{v\ — W(t)) / dv± 


(v±) 


h+l 


'hi-Vool^^ 


(bl - Ko| P + |vi - Wo| 2 + |Ko - W(t)\ P + \Voo - w(t)\ 2 }d Vl 


T 


P+1 


1 .T 


i coCP»(j) + c B C-l,j 

where the last inequality comes from p < 2 and the definition of W. Note that no extra decay can be 
obtained from the vertical component involving v±. Now min (m,q — 1) > — due to (13.61) and (13.71) . so 
that by (13.81) we have a > 1 > 1 — a and 

(p + 1)(1 — a) < pa + 1 — a, for p > 0 

that is, the second term can be absorbed into the first one: 


|Hi(t)| < c G C 


1 


1 


<c G C 


1 


f(i-«)(p+i) 


^( 1 -a )(P+ 1 ) p><r+(l-a) 

which is integrable in t due to (13.51) . 

We now turn our focus to H 2 (t), which is split into two parts H 2 = if 2 i + H 22 as follows. 


(3.10) 


H 21 (t) = dS x L(v\ — W(t)) 

JD(t) JS C 


S ts J{s)ds 


dv 


H 22 (t) = / dS x / L(vi - W(t))/i(T,x(T;t,x, v), v(T;t,x, v))dv 

J D(t) Js c 


where/i(T,x(T;f,x,v), v(T;t,x,v)) = f Q S ts J(s)ds. Now 


|H 21 (t)| < c G C 


0 

r <is. f f 

D(t ) JS C JT 

■ t 


L(v ! - W{t)) 


(s) q ( Vl ) h+1 ( V ± ) l2+1 


dsdv 


(3.11) 




c G C [ ~z f C (| Vl - W(t)\ 2 + | Vl - W(t)\ p ) —ir■ 

JT sy Js c (Vi) 


(VI ) tl+1 (v±) h+1 


dvds 


it c G C / , , n ds it c G C ——- = c G C — 7 — 

Jt (s) T ?- 1 1 ) 


which is integrable over t € (0, 00 ) if (13.61) is satisfied. We note that the constant C in H 2 1 depends on 
Voo. It remains to estimate H 22 , which we put into Lemma 13.31 □ 

Lemma 3.3. Assuming the decay conditions B and iTOI) . we /lave 

1 1 


|H 22 (t)| ^ cgC 


poo 

and so / |# 22 (£)|d£ < oo. 

Jo 


(t) am 

Proof. Because the t ^ 1 case has already been taken care of in Lemma 13.21 we restrict ourselves to t > 1 
in this proof. Recall that we have chosen T = t a < t. For brevity we denote 

x = x(T;t,x,v), v = v(T;t,x,v), 

and 

x* = x(s;T, x, v) = x(s;t,x, v), v* = v(s; T, x, v). = v(s; t, x, v). (3-12) 

By Lemma T3. II with (f, x, v) replaced by (T, x, v), we have 

111 1 


I dd'22 (t) | < cgC 


dS'x 


/£»(t) Js c 2o 


dv / ds L(m — W(i)) 


<s > 9 W) m (m ) il+1 (v ± ) h+1 


< 


c G C f dS x I dv 1 f ds L(vi — W(t))- 

JD(t) Js c Jo 1 


1 1 


1 


iD(t) j s c Jo ( s ) (^i) (m)' 1+1 

We split H 22 into three parts H 22 = H 22 \ + H 222 + H 223 , where H 22 1 is the integral over the set 

{H <i}ns c , 


(3.13) 
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H222 is the integral over 

{|*i| > (t/ 2 )\ Vl - Pool > bT/ 2 } n {h| > 1 }, 

and H223 is the integral over 

il*i| < (t/ 2 )\vi - Fool} n {hi > 1 } n s c . 

We begin with H 2 2i- Solving the characteristic equations (12.131) for x\ with final data at time T using 
m and (13.121) . we have 

\x*\ > |*i — Tv 1 + sriil — c 2 (T — s) > \x\ — Tv\\ — T\v\ — iq| — s|#i| — c 2 (T — s) (3-14) 

where c 2 = 0(cg )■ We estimate two of these terms as follows. On the one hand, considering the time 
interval (T, t) and using the characteristic equation (12.131) , we have 

1 


T\vi - Ui| < CT [ -\qdp< 
Jt \P) 


< C 


(T) q 2 

since q > 2. On the other hand, for t > 1 we have 

\xi-Tv i| > t\Voo -Uil - \xi -X(t) - (T-t)vi \- |X(t) -tV x 

- bT Li, w d "' ir i - i ’.i* 


bT-C -C- 


t 


(T) 9-1 

for some constant C, where we used the fact that 


^ bT-C 


(3.15) 


) dr = C q 


1 


1 


It \(T) 


9-1 


( T ) 


9-1 


dr < C„ 


It (T) 


—j- = C q -~r < (7 ^ c 

q -1 <1 /rpxq-l ^ 1 ^ 


{T )9 


since q > 2 and u(q — 1) > 1. Therefore, combining (13.141) with the two previous inequalities and using 
h | < 1 + c 2 , we have 

|hl > bT~C-c 2 T- h|s + c 2 s > (2c 2 + 2)T-C'-(l + c 2 )T>T-C' 

for 0 < s < T = t a if we choose b > 3c 2 + 2. Plugging the lower bound of |*{| given above into (13.1311 . 
we have 


1-^221 (i) I ^ c gC 


1 (\ Vl ^W(t)\ 2 + \ Vl -W(t)\ p ) 
(s) q 


(T) m (ui )' 1+1 


dsdv 1 


< cgC lrn . m = cgC- 


(T) " {t) am 

Since am > 1, |i/ 22 i(i)| is integrable in t, as desired. 

Now for H 222 , we have, by definition, that |x*| > bT/2. Hence 


1 ^ 222(01 < c G C J J^qds 


1 AW-w{tt + w-wm dv ^ CGC i =CGC _ 1 


vd+l 


(T) r 


(t) c 


(T) m ( Vl y 

which is also integrable in t. 

Finally considering the term if 22 3 and using the characteristic equation (12.131) in the interval (s,T) 
as before, we have 

|a:*| > |*i -Tvi + svi\~c 2 (T-s) > t\Voo - Vi| - |*i — X(t) - (T — t)v i| — \X(t)-tVoo\ -s|hi| -c 2 T + c 2 s 

where again c 2 = 0{cg)- Consideration of the interval (T, t) provides a constant bound of the second 
and third terms on the right side of this inequality. In fact, we already proved in (13.151) that 

\xi - X{t) - (T - t) Vl \ + | X(t) - tV, x \ < C. 

Within the integration region {| \vi — I4o| > |*{|} of H 22 3 , we therefore have 


2 \ V oo ~ i'i| > |au| ^ t\Voo -vi\ — C- c 2 T - s(|f5i| - c 2 ) 
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where C is some fixed constant. Hence 


s > 


1 


\vi\~c 2 1.2 


\V 00 -v 1 \-C-c 2 T}> - 


-c 2 \T-C\> 


c 5 T 


|#i I ^ c 2 


= So, 


1^11 - C 2 

where we have used the definition of S c and where C 5 > 0 provided that we choose b large enough that 
b > 4C + 4c 2 . Moreover, choosing c G small enough, we have by GOD that 


|vi| - c 2 > |ni| - \vi - tii| - c 2 > 1 - Cc G ~c 2 >\. 


Thus from (13.131) we have 


r°° 1 r c 

I# 223 ( 2 )I < c G C / —gds / 

Jso ( s > Jo 


(\vi-W(t)\ 2 + \v 1 -W(t)\ p ) 


(vi) 


/l + l 


dv 1 


€ 


c G C 
Ti- 1 


(*h) 


q -1 


( Vl ) h+1 


(\vi-W(t)\ 2 + \v 1 -W(t)\ p ) dv G 


Thus 


I #223 ( 2 ) | ^ 


c G C 
\q -1 


c G C 


{T ) q ~ 1 (t) a( ' q ~ 1) 

because Zi > q + 1. It is an integrable function of t because a{q — 1) > 1. 

We summarize the estimates of H(t) in the following theorem. 


□ 


Theorem 3.1. Assuming /1 = min (to, g — 1) > we have \H(t)\ ^ c G C(t) a , where a is given by 

(EU). 

Proof. By Lemma 13.21 we have 


|#(2)| < c G C 


< c G C 


1 


(t) 


(l-a)(p+l) 

1 


(2> C 


(t) 


a(q-l) 


1 


(^(l-aXP+i) Mamin(m,g-1) 

The inequality is optimized by equating the two powers. That is, we choose a = (p + l)/(p + 1 + min (to, q — 1)). 
Then we get the stated inequality. If min (m,q — 1) > , we automatically have 


i:niin(m, q — 1) > 1, (1 — a)(p + 1) > 1. 


□ 


4. Effect of the precollisions on the body 

Because we have already estimated the term F 0 q(W) — F 0 (t, X, W ) in Theorem 13.11 we may now focus 
on the estimate of Rw(t ). 

We begin with a couple of observations. First, recall from the earlier discussion in Lemma T2.2I that 
a particle with position x and velocity v at time t may have collided with the body at various earlier 
times t > ti > t 2 > .... Given a particle colliding at time t, we define the infinite collision set Z(t) to 
be comprised of all points (x, v) for which the particle also collides with the body at a sequence of times 
Sj t. If vi(s) denotes the horizontal velocity of such a particle at time s, then 

f W(s)ds= f vi(s)ds 

J Sj JSj 

for the sequence of times Sj —>• t, so that W(t) = v\(t) = v\. Thus such particles comprise a set of 
measure zero and so provide no contribution to the force so that we may neglect them. For each of the 
remaining particles the collision time t is isolated. 

Secondly, away from Z(t), each particle that collides at time t either has a last precollision time 
r(t, x,v) < t or else has no precollision at all in the interval (0, t). Accordingly, we can write the particle 
density at time t as the sum 

f-(t, x, v ) = xi(2, v)/+(t, x(t; t, x, v),v(t; t, x, v)) + xo(2, v)/ o (?)(0; t, x, v)) 


(4.1) 
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where the precollision characteristic functions Xi(L v) and xo(L v ) are defined as follows. Xi(L v) is the 
characteristic function of 

|(£, v) 3(t, x) e (0,£) x D(t) s.t. J W(s)ds = J v x (s;t,x,v)ds, J t;_l (s; t, x, v)| ds < 21 ? | 

with R being the radius of the disk. Moreover, xo(L v ) = 1 — %i(£, v ). The following lemma provides a 
crude estimate on the velocity. 

Lemma 4.1. If x i(t, v) 7 ^ 0; tften 

ill < s "p W,< + ^ 

where (f) s ,t denotes the average of f over the interval [s,t]. 


Proof. Taking the last collision time r = r(f, x, v) before t, we have from (12.1311 
(t - t)vi - c G J J -j^-qdpds < J £i(s;t,x, v)ds < (f - r)ui + c G 


<P)' 


jdpds 


Now 


0,// -j^gdpds ^ C 


ds 


ds 


< 


c Hw^ =c(t - T) ( {s) ~ q+ % 


- ( s ) 9 ” 1 ( i ) 9 - 1 

with C depending on 5 . Hence 

vi(t - t) - c G C(£ - r) <^(s) _9+1 ^ J i>i(s-,t,x, v)ds < i>i(t - r) + c G C(f - r) ^(s} _9+1 ^) 
Since t)i(s; t, x, v)ds = (t — t) (W) t t for xi, we have 

Vi - c G C {(s)~ q+1 ) T t < <W) Tit < v, + c G C {(sy q+1 ) T t ■ 


That is 


Using the fact that 


i»f, (<W%, - C 0 C(< S )-««) T ( ) < »■ « o SUP ( (<W%, + C 0C (W-««) T ( ) . 

sup ((s) _9+1 ) 


c_ 

T,t ^ (t) ’ 


with <7 > 2 , we have 


oc 


inf (W) T t - ^ < vi < sup (W) + 


□ 


We next observe that both the collisions and the flow preserve the product structure of the particle 
density. 

Lemma 4.2. The density f(t,x,v) has the product form /(£, x,v) = a(t,x i,Ti) b(t,x±,v±). In fact, 
b(t,x±,v±) = b o {v±(0;t,x±,v±)). 

Proof. Recall that /o(v) = ao(vi) bo(v±). Next consider any particle (x, v) at a time s for which 
f(s, x, v) = a(s, Xi,vi) b(s, x±, v±). Let t > s be times such that there is no collision in the interval [s, t\. 
Then by (12.1411 we have 

/(s,x,v) = f(s,x(s;t,x 1 ,v 1 ),ir(s-,t,x 1 ,v 1 )) 

= a(s, X!(s; t, xi, Vi), i5i(s; t, x lt V]_)) b(s, x±(s; t, x±,v±),v± (s; t, £_l, v_l))- 

On the other hand, consider a particle that collides at time t and for which the incoming density has 
the form /_(£, x, v) = a(t,x 1 ,^ 1 ) b(t,x±,v±). Then 

/+(£, x, v) = [ K(v 1 -W(t),ui-W{t))f-(t,x,ui,v±)dui 

J (u 1 -W(t))(v 1 -W(t))<0 

= a+(t,xi,vi) b(t,x±,v±) 
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where 

a+(t, X\, Vi) = / K(vi — W(t),ui — W(t)) a-(t,xi,ui) du\. 

J (ui-W(t))(vi-W(t))<0 

Thus the product structure is preserved under both the flow and the collisions. The last statement of 
the lemma is clear because the vertical component is unaffected by the collisions. □ 


Lemma 4.3. As long as 7 and c G are small enough, we have 

\A"f p+1 + 7 + c g) P+1 


(1 + t) p+1 


(4.2) 


\Rw(t)\ < C- 
where C is independent oft,x,c G and A. 

(A 7 p+ 1 +7+cg) P+1 

With a slightly finer proof, one can prove the better rate of decay \Rw(t)\ < C- - ( 1+t y+d -• 

However, by (12.91) we have cr < p + 1 and hence a is a slower decay rate. So we will only bother to prove 
the p + 1 decay rate for \Rw(t)\. 


Proof. Without loss of generality, we may consider only the right side of the body because the left side 
provides the same decay rate with the same proof. We begin by proving an upper bound on /+. We 
claim that 

bo(v±(0',t,x±,v±))a* + < 2C 2 b o (v±(0-,t,x±,v±)) (4.3) 

where 

a* + = sup {a+(r, £; Ui) | £ = X(r), r € [0, 00 ), and m £ [Wo - 3y, + 37 ]} . (4.4) 

Assuming (x, v) ^ Z(t), let to = fo(L x i v ) be the first collision time after time t and r(f,x, v) be the 
last collision time before time t. We use ED to split 


/+(f,x,v) = 


and have 

and 


/ K{v 1 - W(t),ui - W(t)) xi(t,ui,uj_) 

Jui^W(t) 

x/+(t, x(t; t, x, ui,v±), vi (t; t, x, Ui,v±),v±(t-, t, x, m,v±))dui 
+ [ K(v 1 -W(t),u 1 -W(t)) xo(t,U!,v±) 

J Ul ^W(t) 

xfo(vi{0; t, x, ui,vj_), u_l(0; t, x, u 1 ,v±))dui 
IP II 

\vi - Vl(to', t, 27 , Wl)| ^ C G 


(4.5) 


crC CrC 

inf ( w )st 0 - 77T ^ ^i(*o;t,a:i,?;i) < sup (W) st J - 

s ^ t o (to) S^to 


(4.6) 


(4.7) 


(to) ' < to )- 

Here, (14.61) . which we have used many times, comes from (12.131) and (14.71) comes from Lemma [4.11 We 
deduce from (14.61) and (14.71) that 

|vi - W(t)\ < 2(7 + c g ) < 37 (4.8) 

because we assumed 7 > 2 c G . 

We can now estimate (1431) . By Lemma 14.21 

I = [ if(Bl-^((),U 1 -W(t))Xl(t,«I,«l) 

xa+(r, x(t; t, Xi,ui), Vi (r; t, Xi, m)) b Q (vj_(0; t, x±,v±))du! 


and 


II = 


/ K(v\ — W(t),u\ — W(t))\o (t,ui,v±) 

xa o (hi( 0 ; t,xi,ui)) &o(D_l( 0 ; t,x± , v±))dui. 

For J, we use Lemma [4~T1 to deduce 

r W(t) 

K(v 1 — W(t),ui — W(t))du\. 


[\ 

I < b o (v±(0;t,x±,v±))a* + / 

J in: 
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From the definition of W, we deduce 

W(t) < Ko + ^ < V ~ + CA 1 

inf (W),,t > Voo - - f ^ - CA 1 . 

o <s<V t Jo (1 + S) a 

because A^l. Hence 

I < bo(vj_(0; t, xj_,v±))C\A (7 + c G ) 0 + 

because I< is bounded. 

For the second term of dMD, we know as before that |ui(0;t, X(t),u\) — ui| < C. Hence 
II ^ b 0 (v±(0;t,x±,v±))C [ K(v 1 - W(t),u 1 - W(t )) {ui)~ h dux 

since ao(v\) < C(v i) _Zl . Noticing (14.81) . we use condition (12.71) to deduce that 
II < Cb 0 {v±(0;t,xx,v±)) for v x € [V^ - 2j,Voo + 27 ] • 


Thus (14.51) becomes 


/+(f,x,v) ^ b 0 (v±(0-,t,x±,v±)) {C\A (7 + c G ) a* + + C 2 }. 


Using Lemma T4. 21 on the left side of this inequality as well, we have 


a+(r, x\ (t; t, xx,ux), Ui(r; t, x 1 ,u 1 ))b o (y±(0; t, x±, v±)) 
< bo(v±{0;t : xj_,v±)) {("{H (7 + c G ) a+ + C 2 } 


so that 


b 0 (vj_(0;t,x±,vj_))a* + < M+l(0; t, x±, v_x_))CxA (7 + c G ) a* + 

+C 2 b o (v±(0;t,x±,vj_)). 

So choosing 7 so small that A (7 + c G ) < (2Ci) _1 , we have proven the claim (14.31) . 

We now estimate R^{t). By (14.11) . we have 

|/_(f,x,v) - f NB (t,x,v)\ 

= lxi(i; V )/+( T > x(t; x , v), u(t; t, x, v)) + Xo(*, v)/ o (u(0; t, x, v)) - f NB (t,x.,v)\ 

= Xi(*,v) /+ (t, x(t; t, x, v), u(t; t, x, v)) — f NB (t, x, v) 

= Xi(i, v)6 0 (r>j.(0; t, x_l, a+(r, Xi(t; f, Xi, iq), ui(t; t, Xx,vx)) - a 0 (v i(0; t, xi, iq)) 
< Cxi{t,v)b o {v±{0;t,xj_,v±)) 

where in the last line, we used (14.31) . Hence 

|-RwW| < [ [ L(vi-W(t))\f-(t,x,v) - f NB (t,x,v)\dv 

J D(t) Jv,<W(t) 


+ c 
< c 

+ c 


'D(t) Jv x ^W(t) 


L{v r - W(t))xi(t, v)b o (v±(0-,t,xj_,v±))dv 


' D(t) 
rW(t) 

'inf T^t(W) Tt py- 


■W(t) 

inf T ^ t <W> rt -^|y^ 

L(ui — W{t))dv\ 


r r fW(t) 

/ b 0 (v±)dv± L(vi — W(t))dvi 

JD(t) J J 
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Recalling estimate (14.91) . we then have by assumption on K that 

rW(t) 


r w 

\Rw(t)\ < C Im - W(t))\ P dv! 

Jv^-h Pn(s)ds-^P- 


< C 


< c 


Vao-^vWda-^ 

T , 1 f f (Oy + A'] p+1 ) 

00 tj 0 (i + s y 

{Ar/ P+1 + 7 + cg) p+1 


ds — W(t) — 


c G C 


(t) 


p+i 


(l + t)P+ 1 

where in the last line we used the fact that 

1 y 7 + A-/P+ 1 ^ ^1 + A^ +1 

t In 


(s) c 


-ds < C- 


(t) 


□ 


5. Proof of Theorem 12. II 
We begin with the following elementary lemma. 

Lemma 5.1. Suppose 


dY 

dt 


= —b(t)Y + d(t) 


for t > 0, where b(t) >b o > 0 and \d (f)| < Co(l + t) a with a > 1. Then there exists C\ such that 

\Y{t)\ < \Y(0)\e~ bot + C 1 (l + t)~ ,J 

where C\ = 0 (Cq) as Cq —»• 0. 

Proof. Let B(t) = f*b(s)ds. Then ^[e B ^Y (£)] = e B ^d(t) so that 

Y(t) = e- BW Y{0)+e- BW [ e B(s) d(s)ds. 

Jo 

Thus t 

\Y(t)\ < |T(0)|e- bot + [ e - f> o( t - s )( 7 0 (i _|_ s y 7 ds. 

Jo 

Estimating (1 + s) _<T < 1 in [0, i/2] and (1 + s) _,T < (1 + t/2)~ a in [i/2, i], we find that 

\Y(t )| < |W(0)|e- b °* + ^ e ~ bot/2 + ^ (l + 0 

Lemma 5.2. There exists a constant A such that for small enough 7 and cq, we have 

A 7 P + 1 


□ 


\Vw(t)-V 00 \<je~ b ° t + 


(ty 


Y(t) = V w (t) - Voo, b(t ) = 


Foo{W(t)) - FoojVn) 

W(f) - Poo 


In other words, Vw G W. 

Proof. We define 

and 

d(t) = R w {t) - F 00 (W(t)) + F 0 (t,X(t),W(t)). 

We have |T(0)| = 7 by definition. From Theorem 13.11 and Lemma POl we have 

\d(t)\ < c G C(t)~ a + C (Ay p+1 +7 + c G ) P+1 (t)” 1_P • 
Since a < 1 + p by its definition (12.91) , we have 

|d(t)| < C (c G + (A^ p+1 + 7 + c G ) P+1 ) (t)~ a ■ 


Now we apply Lemma Ib.ll to obtain 

\Vw(t) - Pool < 7 e_h °* + C (c G + (^7 P+1 + 7 + c g) P+1 ) (t)~ a ■ 
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We choose A > 2C and c G < \") 2p+1 - We then have 

C (ca + (AY +1 + 7 + c G ) P+1 ) < C 7 P+1 (|7 P + (A 7 P + 1 + \l 2p Y +1 ) < ^7 P+1 
by choosing 7 small enough, since p > 0. That is, 

\Vw{i)-V 00 \< 1 e-^ t + ^-, 

as claimed. □ 


We now have 


dVw 

dt 


E-F 0 (W) 

Voo-W 


(Voo 


Vw) + d(t) 


with Vw{ 0) = Vo and 7 = \Vq — I4o|. Recall that W G W means that W(-) is Lipschitz, W(0) = Vo and 
| W(t) — Wo| < 'ye~ b ° t + ■ Given L > 0, define 


K = {W G W | esssup{\W{t)\ + \W{t)\) < L} 


which is a convex and compact subset of Cb([0, 00 )). Define the operator A on K, by A{W) = Vw- A 
choice of L sufficiently large implies that A : K, —»■ 1C. In fact, 


\v w {t)\ 

Vw(t) 


^ Voo + 37 

« 37 KElv.Xv.M,! f oo(C) + c ((7 + A-y p+1 + C G ) P+I + ca) 


as already shown in the proof of Lemma 15.21 The next lemma will allow us to apply the Schauder fixed 
point theorem to deduce that A has a fixed point, which will complete the proof of existence in Theorem 

m 


Lemma 5.3. If we provide /C with the topology of Cb([0, 00 )), then A is a continuous operator. 


Proof. We follow the proof in (9J Lemma 7.1], slightly modified. The term H = Hw = Fo o(W) — 
Fo(t,X,W) involves no collisions at all, while the recollision force satisfies |i?w| ^ G 7 (t)~ a ^ C 7 . 

For any t > 0, we define B^(t) as the set of (x, v) such that the trajectory passing through (f, x, v) 
has collided with the body at least N + 1 times in [0,f], and we define A§y(t) as its complement. Using 
the estimate (14.211 that sup te r 0 oo ^ vvew |Rw(0l ^ G 7 and iterating this estimate N times, we have 

sup \R w (t.;Bw(t))\ < (C'y) 1 *. (5.1) 

te[0,oo), w ew 


For any e > 0, this expression is at most | by choosing 7 < 1/C and N = N e large enough. 
Now let Wj —> W in K.. Given any time S > 0, we may write 


sup | R Wj {t)~Rw(t) 

0<£<oo 


< 


sup \R Wj {t) - R w {t)\+ sup I R w {t\Bw{t)) 
S<t< 00 0 <t<5' 


sup 

0<t<5 




+ sup 
0 <t<S 


Rwj (t; Aw s (t)) — Rw{t’i Aw(t)) 


= I + 11 +III + IV. 


By estimate (USD, we may choose S = S e so large that |/| < e/4. By estimate (15.111 . we have \II\ + \III\ < 
2e/4. 

Now in IV there are no more than N collisions. Therefore we can express both terms in TV as iterates 
of at most N integrals by repeated use of the collision boundary condition. The resulting finite number 
of iterated integrals contain Wj in a finite number of places in the expression Rwj (t; A§y. (t))■ (See 
Lemma [2.21 1 Therefore they converge as j —» 00 to the same expression with Wj replaced by W, and 
the convergence is uniform for t G [0, S']. Thus we can choose j so large that \IV\ < e/4. So we conclude 
that Rwj{t) —>• Rw(t ) in Cb([0, 00 )). 

It is also clear that 

Hw, -> Hw 
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in /C because no collisions are involved. It follows from the equation (12.181) for Vw(t ) that A is continuous 
in the topology of Cb([0, oo)). □ 

So far we have proven the existence part of Theorem 12.11 We now prove that every solution in this 
sense satisfies (12.101) with a simple argument. In fact, consider any solution 

no) - v c 


(m/) 


that satisfies 


Then 


V(t) - V« 


= 7 < 7 + Aj p+1 . 

,-M ^ A ^ P+1 


< 7e 


(ty 


(5.2) 


for small t. Now suppose we have equality at some later time. Let T be the earliest such time. Then 

(v{t), /) is a fixed point of the mapping 

□ 


the existence of such a time T contradicts Lemma 15.21 because 
A. Thus (15.21) is valid for all t < oo. 
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